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Abstract. The thermoelectric effect in a quantum dot (QD) attached to two leads 
in the presence of microwave fields is studied by using the Keldysh nonequilibrium 
Green function technique. When the microwave is applied only on the QD and in the 
linear-response regime, the main peaks in the thermoelectric figure of merit and the 
thermopower are found to decrease, with the emergence of a set of photon-induced 
peaks. Under this condition the microwave field can not generate heat current or 
electrical bias voltage. Surprisingly, when the microwave field is applied only to one 
(bright) lead and not to the other (dark) lead or the QD, heat flows mostly from the 
dark to the bright lead, almost irrespectively to the direction of the thermal gradient. 
We attribute this effect to microwave-induced opening of additional transport channels 
below the Fermi energy. The microwave field can change both the magnitude and the 
sign of the electrical bias voltage induced by the temperature gradient. 
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1. Introduction 

Owing to the possibility of high heat-voltage conversion efficiency, thermoelectric 
phenomenon in nanoscale solid-state materials has become an active research area 
in recent years[l,2]. In the linear response regime, i.e., the temperature gradient in 
the device AT approaches to zero, the thermoelectric efficiency is measured by the 
dimensionless thermoelectrical figure of merit ZT, and is usually smaller than one in 
bulk materials [3]. This hinders the wide usage of thermoelectric effect in commercial 
applications. The reason behind the small ZT in bulk materials is the Wiedemann- 
Franz law[4], which can be violated in nanoscale materials, for instance due to the 
Coulomb blockade effect [5]. In the nonlinear response regime, i.e., finite AT case, 
the thermoelectric energy conversion efficiency is characterized by the bias voltage 
AV generated by AT [6]. In both of linear and nonlinear cases, the thermoelectric 
efficiency in nanostructures may be much higher than that in bulk materials because of 
the reduced dimensionality [7-9] , and more interestingly, depends on the quantized energy 
levels that can be modulated by gate voltages. Since the breakthroughs in experimental 
work [10- 14], large figure of merit exceeding one or higher and large AV were frequently 
reported in superlattices, quantum wires, quantum dots, and carbon nanotubes, etc. 
Theoretically, impacts of Coulomb blockade[5-9, 15-20], Kondo[21-24], and Fano[25,26] 
effects on the thermoelectric properties were extensively investigated. Applications 
such as nanoscale refrigeration [27], thermal rectifier [18, 28, 29], thermal transistors [6], 
and thermal memory and logic gates [30] were proposed. 

The above-mentioned works mainly focused on the thermoelectric efficiency under 
time-independent fields. Very recently, heat transport with photons has also been 
studied. It was demonstrated that heat can be conducted by photon radiation at very 
low temperature when the phonons are frozen out [31]. Quantum circuits sandwiched 
between two reservoirs and threaded by electromagnetic fluctuations (photons) were 
proposed to operate as heat transistors or electron coolers [32, 33]. In a double-quantum- 
well structure, it was shown that hotter electrons in one reservoir may be replaced 
by the cooler ones in the other reservoir with the help of photon by adjusting the 
wells' energy levels[34]. Electron heating in two-dimensional electron system has 
also been theoretically and experimentally investigated [35]. In addition, Lau et al. 
have calculated the properties of coherent radiative thermal conductance in photonic 
crystals[36]. Moreover, enhanced thermopower (Seebeck coefficient) induced by a time- 
dependent gate voltage was proposed in a device of a quantum dot (QD) coupled to two 
metal leads [37]. 

As a zero-dimensional system, QD is a promising thermoelectrical material and was 
intensively studied in many of the above mentioned works. Until now, photon-assisted 
heat transport in such a device has seldom been addressed, although it has been studied 
in the context of electron transport through QDs long time ago[38,39]. In this paper, we 
study the photon-mediated thermoelectric effect in a QD attached to two leads in both 
linear (AT — > 0) and nonlinear (finite AT) regimes. Compared with previous studies on 
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quantum circuits [32, 33], we deal with a relatively high temperature situation. Our main 
results show that when the microwave field is applied only on the QD, the magnitude of 
the figure of merit is suppressed because of the photon-induced sub-channels in the dot. 
When the microwave is applied only on one lead, large heat current and electric bias 
are generated. We find that the heat current generally flows from the dark lead into the 
bright one, except for the situation in which the dot level is around the electron-hole 
symmetry point. In the present paper, the contributions from phonons are neglected 
for the sake of simplicity. 

2. Model and method 

The present system can be described by the following Hamiltonian[39,40]: 

+ ^(W^A + H.c), (1) 

where c\ j3a (c^) is the creation (annihilation) operator of the electrons with momentum 
k, spin a and energy Skpif) in the j3 (= L, R) lead; d\ {d a ) creates (annihilates) an 
electron with spin a, energy Edit) in the QD. U denotes the intradot Coulomb interaction. 
The last term in Eq. (1) represents the dot-lead coupling with coupling strength Vpa. 
The time-varying energy spectra originate from the applied microwave fields and take 
the forms of Ekp(t) = E^p + Apcosut and Edit) = Ed + A^coswt, where £ k p and Ed 
are the time-independent single electron energies without microwave field, and Ap/d, 
uj are respectively the microwave strength and frequency applied on the leads or the 
dot [39,40]. 

Time- dependent electric and heat currents from lead (3 to the QD can be calculated 
from the evolution of the total number operator of the electrons in the lead (in unit of 
h=l)[39], 

( Qui) ) = ( Z(t) - E F {t) ) C ^ Cfc/3CT) ' (2) 

where Ep(t) = e[/Xi(t) — is the Fermi level. The total electric and heat currents 

are J(t) = Y.A J La{t) - Jna{t)} and Q(t) = J^AQ^) ~ QrA 1 )], respectively. In this 
paper we are interested in the time averaged electric and heat currents, which can be 
derived following the standard Keldysh Green function technique as [39]: 

(o) = h^h^ I d£ ( e- e f ) 

x [f L (e)lm(A La (e,t)) - f R (e)lm(A^(e,t))}, (3) 

where the line-width function = ^Ylk^H^d^Lfrnd^i 6 - ~ Ek ) * s assume d to 

be independent of the energy under wide bandwidth approximation. //3(e) = [1 + 
e ( e- M/3)/fcB T /3]-i j s the Fermi distribution function of lead (5 with chemical potential /ig, 
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temperature Tp and Boltzmann constant k B . The quantity Ap a (e,t) is defined through 
the Green function as 

A Pa (e,t)= [ dtiGlfatJeM-Hh-Q-i [ A (r)dr], (4) 

J — oo J t 

where G r a (t,ti) = —i9(t — t')({d a (t) : (£')}) is the retarded Green function, which can 
be obtained by the equation of motion technique following the processes in Ref . [40] : 

G r a (t,t')= -ifl(t _*')){(! _„ ff )expH j\ d {T)dr-^{l-n^){t-t')) 

f* r 

+ n^exp(-ij [e d (r) + U]dr - -n„{t - t ))}. (5) 

It should be noted that the adopted decoupling approximation to higher-order 
many-particle Green functions is sufficient for electronic transport in the Coulomb 
blockade regime, but not for higher-order tunneling process such as the subtle Kondo 
effect, which is beyond the consideration of this paper. In the above equation, 
r = T L + T R . Substituting the expression of the retarded Green function into Eq. 
(4) and carrying out the integrations, A^ a (e,t) becomes 

,A d - A? A/3 - A c 
)Jk'\ 



k,k' 



x -f ^ n<J + n<T 1 (6) 

e — Ed — k'u + iT(l — rba)/2 e - £ d - U - k'u + iTn- a j1 h 1 ' 

where are Bessel functions of the first kind. The occupation number n 5 in the above 
equations needs to be calculated self-consistently from the equation 

n- c = QmG<(t,t)) = J ^^2Me)mA^(e,t)\ 2 ), (7) 

where G<(t,t) is the lesser Green function that can be calculated by the Keldysh 
equation with the help of the retarded Green function. Once Ap a {e,t) is determined, 
the currents can be calculated from its time-averaged form: 

(A ^t)) = Y J Jl{^ 1 ) 
— ui 

k 

x r \ ~ n ° | ^£ 1 /g\ 

[ E-E d -kuo + iT(l -n 9 )/2 e - E d -U - koo + iTn^/2 1 ' 1 ; 

3. Results 



Let us start by considering the case where the microwave field is applied only on the 
dot, i.e. A L — A R = and A d ^ 0. In the following numerical calculations, we set 
h — 1 and choose the microwave field frequency u = 1 as the energy unit. In the linear 
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response regime (T L — T R — T) the electric current in Eq. (3) can be written in the 
following compact form: 

J = - 2 -T f de[f L (e)-f R (e)]T(e), (9) 



h „ 

where the transmission coefficient T{e) is 

i- L + i- R , U 

X K-Ed-ku + iYil-n^/2 + e - e d - U - ku + iTn & /2^' ^ 

We subsequently introduce the integrals I n {T) with n = 0,1,2, I n (T) = 
— {2/h) f e n (df /de)T(e)de. The linear conductance G, thermopower S, and the thermal 
conductance n can be expressed as 

G = e 2 / (T), 

S = -h(T)/[eTI (T)\, 

k =(l/T)[/ 2 (T)-/ 1 2 (T)//o(T)]. (11) 

Neglecting the thermal conductance from the phonons as is assumed here, the figure of 
merit ZT is given by 

ZT = GS 2 T/k. (12) 

Figure 1 shows the behaviors of G, k, S, and ZT versus dot level Ed for different 
microwave field strength A d . For A d = 0, the linear conductance G in Fig. 1(a) shows 
two peaks respectively positioned at e d = and —U, exhibiting the typical Coulomb 
blockade effect. As for the thermal conductance k in Fig. 1(b), two smooth shoulders 
emerge around the two main peaks. When the microwave field is applied on the dot, 
photon-induced peaks emerge in G at ± ku and —U ± kcu, where k is an integer. 
Meanwhile, the height of the main peaks is reduced. The reason is that the microwave 
field excites more conduction channels in the dot, and then the transmission probabilities 
through the channels of e d = and —U are shared and suppressed. As for the thermal 
conductance k, the height of the main peaks first increases and then decreases with 
the increase of the microwave field strength. This is because k is determined by both 
the transmission probability and the heat transferred by each electron. Contrary to the 
situation in G, each main peak in k around e d = or —U is contributed from more 
than one conduction channels, which can be seen by comparing Fig. 1(b) with 1(a). 
Thus, their height may be enhanced even if the average tunnelling probability of each 
channel decreases. Here the function of the microwave field is analogous to the changing 
of temperature in the leads[18-20], which adjusts the numbers of the electron and hole 
around the Fermi level to participate in transport. With increasing temperature T, 
the average tunneling probability will be monotonically suppressed but the heat carried 
by each electron increases. As a result, the height of the peaks in k will exhibit a 
nonmonotonic dependence on T [19, 20], which is different from the behavior of the 
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Figure 1. (Color online) Electrical conductance G, thermal conductance k, thermal 
power S, and figure of merit ZT as functions of dot level £d for different microwave 
field strength with Al = = 0. Other parameters are U — 3w, ksT = O.loj, and 
= Tr = 0.05w. The presence of a microwave field is seen to decrease the strength 
of the main peaks of the thermoelectric figure of merit, and a set of additional peaks 
in multiples of the microwave frequency emerges. 



peaks in G. We also observe that a peak in the thermal conductance develops at the 
electron-hole symmetry point —U/2 due to the application of the microwave field. This 
may be attributed to the photon-induced conduction channels around this point, which 
can be seen from Fig. 1(a). 

The thermopower S in Fig. 1(c) has three zero points individually at Ed = —U, 
—U/2 (electron-hole symmetry point), and 0. The magnitude of S oscillates between 
these three points with positive and negative values, which can be understood as 
follows [19, 20]: In the presence of temperature difference that induces the thermoelectric 
phenomenon, there are more electrons above the chemical potential fi = in the hotter 
lead, and at the same time more holes in the cooler lead below \x. When the dot energy 
levels Ed or Ed + U is above /i, more electrons will transport from the left lead of higher 
temperature to the right lead of lower temperature than the holes tunnelling from the 
opposite direction, inducing a positive voltage drop AV and consequently a negative 
thermopower S = —AV/AT (S is negative in the figure as in unit of fce/e). When the 
energy levels are aligned to the chemical potential or the electron- hole symmetry point, 
the voltage drop induced by the electrons are cancelled by that of the holes. Thus the 
thermopower is zero. When Ed or Ed + U is below the chemical potential, more electrons 
tunnel from right to left and S changes sign. The above explanation also holds true when 
the microwave field is applied. Now a set of sub-channels are opened by the photons, 
inducing more zero points in S. Consequently, the applied microwave field can change 



Microwave-mediated heat transport in a quantum dot attached to leads 



7 



the sign of S, which may be used to detect the properties of the microwave field. Fig. 
1(d) shows that the magnitude of the main peaks in ZT is reduced by the microwave 
field with the emergence of a set of shoulder peaks. This is because of the combined 
effect of reduced linear conductance and enhanced thermal conductance. Although the 
original peaks height is lowered by the microwave field, there are several new peaks that 
can reach quite large value, and may be useful for energy conversion. 

Figure 2 shows the temperature-dependence of ZT for different microwave field 
strength A d . To acquire a large ZT value, the dot level is fixed at -0.25a;, where 
the main peak in ZT is centered. It is shown that the ZT value first increases and 
then decreases after reaching a maximum value with the increase of the temperature, 
which is consistent with previous result [18-20]. The magnitude of ZT is determined 
by the combined behavior of G, k, and S. Detailed dependence of these quantities on 
the temperature can be found in Refs. 19 and 20. For ksT <C U, the ZT value is 
considerably enhanced due to the Coulomb blockade effect, and optimal thermoelectric 
efficiency can be achieved in this region. With the increase of the temperature, the peaks 
in ZT is suppressed as the influence of Coulomb blockade effect on electron transport 
becomes less important. In the presence of the microwave field, the line-shape of ZT 
is reserved but with reduced strength. This is consistent with the above result. With 
the increase of the microwave field strength, the maximum ZT value shifts to lower 
temperature regime. 




Figure 2. ZT as a function of ksT for different microwave field strength = 
0, cj/2, uj, 3iv/2 with A^ = A^ = 0. Other parameters are as in Fig. 1. The effect of 
microwave field applied only to the quantum dot is a string reduction of ZT. 

When the microwave field is asymmetrically applied to the structure, heat current 
and bias voltage will be generated under the condition of zero electrical current (this 



Microwave-mediated heat transport in a quantum dot attached to leads 



8 



condition is easily met in thermopower experiments [1,2]). From Eqs. (3) and (8), now 
the heat current is mainly determined by the difference between the microwave field 
strength in the dot A^ and that in the leads £±l/r- We thus consider the case of a 
field applied only on the leads. For comparison, we first present the results without any 
microwave in Fig. 3(a-b). It shows that the heat current always flows from the hotter 
lead to the cooler one under a thermal bias AT. It has a main resonance when the dot 
level is located on the electron- hole symmetry point, and two dips when e d and Ed + U 
are aligned to the chemical potential /j, — 0. The reason behind this is as follows: When 
Ed and Ed + U are aligned to the chemical potential, the energy difference between the 
electrons from different leads, which are tunnelling in opposite directions through these 
levels, is relatively small, resulting in small heat current. But when the dot level is 
at the symmetry point Ed = —U/2, electrons in the left lead tunnel to the right lead 
through the level of Ed + U = U /2, which is above the chemical potential and have larger 
thermal energy. Meanwhile, equal number electrons tunnel from the right lead through 
the level of —U/2 to the left one with less thermal energy. In other words, holes tunnel 
from the left lead to the right lead through —U/2. The contributions from the electrons 
and holes to the heat current add constructively, giving rise a high peak. Other smaller 
resonances in Fig. 3(a) can be explained by similar reasons. The bias voltage AV in 
Fig. 3(b) has three zero points at energy levels of Ed = 0, — U and —U/2, because 
where equal numbers are flowing from opposite directions. Apart from these levels, it 
oscillates between negative and positive values. The reason is just the same as that of the 
thermopower in Fig. 1(c). If there is microwave field applied on the QD as shown in Fig. 
3(c-d), the trends of the heat current and electrical bias voltage are preserved, except 
for the emergence of a set of smaller peaks. The heat current magnitude is enhanced 
because there are more carriers participating in transport through the photon-induced 
sub-channels. 

Let us now consider the situation when the microwave field is applied only to the left 
lead. Fig. 4(a-b) shows the heat current and the induced electric voltage as functions 
of the dot level for various temperature differences with A^ = A# = and = w. 
Compared to the cases in Figs. 3(a-b) and 3(c-d), the magnitude of the heat current 
is remarkably enhanced as the dot levels Ed and Ed + U are slightly higher and lower 
than zero. For positive thermal bias ksAT > o, the heat current can flow from the 
right cooler (dark) lead to the left hotter (bright) lead except for £ d is in the middle 
of the Coulomb blockade regime, which is shown in Fig. 4(a). For negative thermal 
bias k B AT < 0, the heat current always flow from the right hotter (dark) lead to the 
left cooler (bright) one (see the solid and dashed lines in Fig. 4(a)). As indicated by 
the behavior of the electric bias AV^ in Fig. 4(b), electron transport behavior under 
both positive and negative thermal biases resembles that of positive thermal bias case 
that shown in Fig. 3. Meanwhile, the magnitude of the electric bias AV^ of positive 
thermal bias is enhanced because of the photon-electrical effect. When the microwave 
field is applied on the left hotter lead, the photon-induced peaks in the bias voltage are 
smoothed off due to thermal effects. The above trends of the heat current and AV hold 
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Figure 3. (a) Heat current Q and (b) bias voltage AV induced by different thermal 
bias AT as functions of dot level Ed in the absence of microwave field. The temperatures 
of the left and right leads are T L =T + AT/2 and T R = T- AT/2 with k B T = 0.25a;. 
Other parameters are as in Fig. 1. (c-d) same as (a-b) for different thermal bias and 
in the presence of various microwave strength A^ (detailed in the legend), keeping 
Al = Ah = 0. Both the heat current and the induced voltage exhibit a decrease 
of their maximum values, and the appearance of additional microwave-induced peaks. 
The heat always flows from the hot to the cold lead. 



true even for zero thermal bias case. 

We now qualitatively discuss the impact of phonon on the figure of merit and the 
heat current. At low temperature regime, we can add to the thermal conductance k in 
Eq. (12) a phonon contribution which is K p h = ?>TT 2 k 2 B T /h [16]. Then the magnitude of 
ZT will be reduced correspondingly. But the trend of ZT can be unchanged [19]. In the 
non-linear case, phonon peaks in the heat current may also emerge, which should be an 
interesting topic. Nevertheless, recent studies indicated that the heat current delivered 
by phonons can be blocked by particular device design [19], and moreover, K p h can be 
small depending on the phonon spectrum [41,42]. 

4. Discussion and summary 

The result presented in the last section is rather surprising: although energy (via the 
microwave field) is pumped into the bright lead, still heat flows from the dark to the 
bright lead. The explanation for this effect we point out that the microwave field shifts 
the spectral weight of the electrons in multiples of the photon frequency, thus effectively 
reduces the density of states (DOS) in some places and increases it in others. In this case, 
the effective reduction of the DOS is most pronounced near the particle-hole symmetry 
point, and the weight is shifted below it. The result is electrons, which are flowing from 
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Figure 4. Heat current Q (a) and bias voltage AV (b) as functions of dot level 
£d for different thermal bias ksAT = — 0.2w, — O.lw, 0, O.lw, 0.2w. A microwave field 
Al = u> is applied only to the left lead. Other parameters are as in Fig. 3. Due to the 
asymmetry in the microwave field, heat current flows from the dark to the bright lead 
almost regardless of the temperature gradient, except at the vicinity of particle-hole 
symmetric point Sd = —U/2. 



the bright to the dark lead, carry negative energies, because they are placed below the 
Fermi energy. 

To understand this unique situation it is useful to turn to the analogy with the 
non-interacting dot where U — 0. In that case, the dot is a simple resonant level, and 
consequently Im^A^e, i)) is simplified, 

M^M))^(^) ( ^_ r M2 + r2 . (13) 
Substituting it back to Eq. (3) one obtains for the heat current 

Q = l^S £d£ {s-e d -kuy + T* ti^ J ^ A ^ ~ f^Vm] -(14) 

The Bessel functions Jfc(x) are decaying functions of k for x < 1.434. Since the 
Fermi functions decays above the Fermi level, the main contributions to the integral 
above are k = and k = — 1, resulting in 
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« ~ I J £d£ {£ -^)2 + T 2 [h(e)J 2 (A L /u) - f R (e)] 

+ ll £de l I », f2 /l(^(AlM 

h J (e - s d + u) 2 + T 2 

= ^J e/(e)cfe. (15) 

In the absence of the microwave field, the transmission is centered around Ed and 
has a single peak, but additional conduction channels open when the microwave field 
is applied. This can be seen in Fig. 5 where 1(e) (solid line), which is proportional to 
the transmission function, and the heat current integrand e/(e) (dashed line) plotted 
as functions of the integration variable e, where we have chosen the dot level Ed = 
for clarity. 1(e) shows an additional peak at negative energies, which induces electron 
current from the bright to dark leads. The heat current is thus in the opposite direction. 
Note also that the weight of the negative peak is larger than that of the peak near Ed, 
resulting in a net negative heat current. 

02 
0.1 

m 
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-0.2 

-2.0 -1.5 -1.0 -0.5 0.0 0.5 

e 

Figure 5. (Color online) 1(e) (solid line), which is proportional to the transmission 
function, and the heat current integrand sl(e) (dashed line) plotted as functions of 
the integration variable e for Ed = 0. An additional peak appears in /(e), and when 
weighted with energy el(s) results in negative heat current, i.e. from dark to bright 
lead. 

In summary, microwave-assisted heat transport in a quantum dot is investigated in 
the framework of nonequilibrium Green's function method. It was found that the figure 
of merit can be suppressed by a microwave field applied on the dot. In the completely 
asymmetric case, i.e., the microwave field is applied only on one of the leads, the heat 
current generally flows from the dark lead to the bright lead, weakly depends on the 
direction of the thermal bias, and can be reversed only when the dot level is close to 
the electron- hole symmetry point. The results we have presented here can be verified 
in experiments similar to those of Ref. [38] in which a temperature difference is applied 
between the two leads as well as a microwave field and gate voltage. We note that 
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here we only deal with the fixed microwave frequency case at steady state. One can 
expect frequency-and time-dependent heat transport will bring about more interesting 
phenomena, which are beyond the scope of the present paper, and research in that 
direction is currently underway. 
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